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Abstract. Using the methods of quantum trajectories we study numerically the phenomenon of quantum 
synchronization in a quantum dissipative system with periodic driving. Our results show that at small val- 
ues of Planck constant h the classical devil's staircase remains robust with respect to quantum fluctuations 
while at large h values synchronization plateaus are destroyed. Quantum synchronization in our model has 
close similarities with Shapiro steps in Josephson junctions and it can be also realized in experiments with 
cold atoms. 

PACS. 05.45.Xt Synchronization; coupled oscillators - 03.65.Yz Decoherence; open systems; quantum 
statistical methods - 74.50.+r Tunneling phenomena; point contacts, weak links, Josephson effects 



Since 1665, when Christiaan Huygens discovered the 
synchronization of two maritime pendulum clocks pQ (see 
■ for historical survey and modern experiments), it 
has been established that this universal nonlinear phe- 
nomenon appears in an abundant variety of systems rang- 
ing from clocks to fireflies, cardiac pacemakers, lasers and 
Josephson junction (J J) arrays 0. Various mathematical 
tools have been developed to analyze synchronization in 
simple dissipative nonlinear models with periodic driving 
and complex ensembles of nonlinear coupled oscillators. 
Such pure mathematical concepts as Arnold tongues in 
the circle map £y found their experimental implementa- 
tions with Shapiro steps |B and JJs synchronization [5]. 
However, till recently the treatment of synchronization has 
been mainly based on classical mechanics even if JJs have 
purely quantum origin 

A significant progress in the theory of dissipative quan- 
tum mechanics has been done in Ref.|7|- It was further 
developed by various groups as reviewed in [H] . Nowadays 
this research line is getting more and more importance 
since technology goes on smaller and smaller scales where 
an interplay of dissipative and quantum effects becomes 
dominant. A typical example is given by small size JJs. 
Here dissipative effects are always present even if in cer- 
tain cases skillful manipulations allow to realize long term 
coherent Rabi oscillations [§]. 

These reasons led to a significant number of analyt- 
ical studies of quantum tunneling in presence of dissipa- 
tion (see e.g. |7H51l 11)111 11112) ). However, analytical methods 
have serious restrictions in a strongly nonlinear regime 
typical of synchronization Due to that in this pa- 
per we perform extensive numerical studies of quantum 
synchronization following precisely a transition from clas- 
sical to quantum behavior changing effective dimension- 



less Planck constant h by three orders of magnitude. The 
quantum dissipative evolution is described by the master 
equation for the density matrix p written in the Lindblad 
form 13 . To perform simulations with a large number of 
states TV cx 1/H we use the method of quantum trajectories 
|14U15j which allows to reduce significantly the number of 
equations compared to direct solution of the master equa- 
tion (N instead of TV 2 ). In this approach one quantum 
trajectory can be viewed as one individual realization of 
experimental run |16j . 

For our studies we choose a model which in the classical 
limit is closely related to the circle map |J which gives a 
generic description of Arnold tongues and synchronization 
of system dynamics with an external periodic driving 3 . 
In absence of dissipation the evolution is described by the 
Hamiltonian of kicked particle falling in a static field /: 

H =p 2 /2 - fx + Kcos x S(t-m), (1) 

m— — oo 

with the usual operators x and p = hn = —ihd/dx. At / = 
this Hamiltonian corresponds to the kicked rotator ^7] , 
a paradigmatic model in the fields of nonlinear dynamics 
and quantum chaos. At non-zero / the system has been 
built up in experiments with cold atoms placed in laser 
fields and falling in a gravitational field ^Hj- Nontrivial 
quantum effects of the static force are analyzed in |19| . 
Here we consider the dynamics of the model in presence 
of additional friction force F = —g 2 p. Experimentally such 
a force can be realized by the Doppler cooling 1201 ■ The 
model can be also implemented with JJs as described in 
PT] . In this case effective kicks are created by an external 
ac-current source, / andp are proportional respectively to 
dc-current and voltage drop across JJ, while the friction 
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force F naturally appears due to finite circuit resistance. 
The expression of K and ft via JJ parameters is given in 

eh 

The classical dynamics with friction can be exactly 
integrated between kicks that gives a dissipative map 



P = (1 - l)p + (1 - 7)^sin x + f-f/g 2 , 

x = x + jp/g 2 + (jK/g 2 ) sin x + f{g 2 - ~/)/g 4 , 



(2) 



where bars note new values of variables after one map it- 
eration and 1 — 7 = cxp(— g 2 ). Up to parameter rescaling 
and shifts in x,p, produced by static force, the map 
has the form of Zaslavsky map |221 • Due to contraction in 
p the dynamics in phase variable x is close to the circle 
map x — x + K e f f sin x + v 00] and demonstrates devil's 
staircase structure in the dependence of average momen- 
tum P on / (Fig. ^ top). Steps near rational rotation 
numbers P/2-k correspond to synchronization with exter- 
nal periodic driving inside Arnold tongues. In average the 
momentum P = f/g 2 as it should be in an equilibrium 
between the external and friction forces. 

The corresponding quantum dissipative dynamics is 
described by the master equation in the Lindbald form 

1X3 : 



(3) 



where p is the density operator, { , } denotes the anticom- 
mutator, are the Lindblad operators, which model the 
effects of the environment. Following [23] we assume the 
Lindblad operators in the form (p — 1,2): 



Li = .9E„v / ^" 
^2 = g E n 



1 In) (n 
n) 



1 



1|, 
-n 



(4) 



These operators act on the bases of 27r-periodic eigenstates 
of operator n and in the regime of weak coupling and 
Markov approximations describe the dissipation force F = 
—g 2 p induced by a bosonic bath at zero temperature. As 
in |2.'{j the numerical simulations of quantum jumps are 
done for one quantum trajectory using the so-called Monte 
Carlo wave function approach |lfi| . The additional term 
with the constant force / is exactly integrated between 
jumps leading to a drift of wave function amplitudes in 
the space of momentum eigenstates n. The total number 
of states TV is fixed by a condition of keeping all states 
with probabilities higher than 10~ 7 . 

The numerical results for the quantum devil's stair- 
case are shown in Fig. ^ f° r various ft values The 
dependence is similar to those of I — V characteristics 
shown in [5] with P cx V, / cx /. At small ft values 
(ft = 0.012) the steps remains stable with respect to quan- 
tum fluctuations. As a result, the rotation frequency is un- 
changed when / is varied in some finite interval. This can 
be viewed as the manifestation of quantum synchroniza- 
tion which takes place inside quantum Arnold tongues. 
For larger ft = 0.05 small steps start to disappear and at 
relatively large h = 0.5 the dependence f — P becomes 
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Fig. 1. Dependence of the average momentum P on static force 
/ at K = 0.8 for 7 = 0.25 (left column) and 7 = 0.05 (right 
column); P is computed over t = 500 map iterations. From top 
to bottom: classical case at h = 0, h = 0.012, ft = 0.05, ft = 
0.5. Initial conditions correspond to one classical trajectory at 
x — 0, p/2-K — 0.38 for classical dynamics. For the quantum 
evolution one quantum trajectory is taken at the same x,p 
position with the wave function in the form of minimal coherent 
state at given ft. 



smooth so that a quantum particle slides smoothly un- 
der static force /. It is interesting to compare the cases 
of strong 7 = 0.25 and weak 7 = 0.05 dissipation (Fig. ^ 
left and right columns). At strong 7 the contraction in p 
is rather rapid and the behavior is similar to the case of 
the circle map U and P grows monotonously with /. At 
weak 7 there are many different attractors in the phase 
space and a classical trajectory may jump between them 
rather irregularly with variation of / or initial conditions. 
At small ft a quantum trajectory reproduces this behav- 
ior of steps overlap rather accurately. We note that very 
similar structure of steps overlap is clearly seen in exper- 
imental data shown in Fig. 2 of [o]- 

To understand the properties of dissipative quantum 
dynamics we analyze the structure of wave functions as- 
sociated with quantum trajectories of FigQ] It is known 
that dissipation in quantum systems leads to a collapse 
of wave packet 25 ,26 27 . In agreement with these results 
and recent studies of Zaslavsky map [22! we find that the 
wave function in our model collapses on a compact packet 



O. V. Zhirov, D. L. Shepelyansky: Quantum synchronization 



3 




1 



4 0.0 0.2 



0.4 



0.6 . O.i 



Fig. 



f f 

2. Dependence of rescaled dispersion o^/ft 1 / 2 of wave 



packet in x on force / for the cases of Fig.l at 7 = 0.25 (left 
panel) and 7 = 0.05 (right panel). Symbols show data at ft — 
0.012 (black points), ft = 0.05 (open circles), ft = 0.5 (+). 
Dispersion a x is averaged over t = 500 map iterations for one 
quantum trajectory with initial state of Fig. Q 



of width a x , a p (dispersion) in coordinate and momen- 
tum respectively (e.g. \i/>(x)\ ~ exp(— x 2 /(4a 2 ))) The 
width a x depends nontrivially on / as it is shown in Fig [2 
Inside the steps with synchronization the value of a x is 
strongly reduced while in the sliding regime between steps 
its value is significantly enhanced. In the limit of small ft 
where quantum synchronization takes place the rescaled 
width a x /Vh shows approximately the same functional 
dependence on / for all ft. This functional dependence on 
/ is completely changed at large ft where synchronization 
is absent and where a x is not very sensitive to /. For small 
7 (FigEI right) the dependence on / becomes more irreg- 
ular in agreement with more complicated step structure 
seen in Fig^ (right). We note that at the same time a p 
remains practically insensitive to variation of /. This is 
related to the fact that in momentum p the dynamics is 
close to a simple contraction while in coordinate x the 
nonlinear dynamics significantly depends on the system 
parameters (see discussion below). 
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Fig. 3. Dependence of dispersion a in x (a x ) and p (<j p ) on 
ft for K = 0.8 at 7 = 0.25 (left) and 7 = 0.05 (right); a is 
averaged over t = 5000 map iterations. Left panel: / = 1.1 
(points for a x and x for cr p ), f = 1.9 (open circles for a x and 
+ for a p ). Right panel: / = 0.21 (points for a x and x for a p ), 
f = 0.37 (open circles for a x and + for a v ). Dashed line shows 
the dependence a = (h/2) 1 ^ 2 . Initial state is as in Fig. 



in the phase space and another is taken on a slope when 
the classical dynamics has attractor in a form of invariant 
curve as it is shown in Fig.0]for 7 = 0.25. Clearly a x is sig- 
nificantly larger in the case of invariant curve than in the 
case of fixed point. Intuitively it is rather natural since in 
the first case variation of phase x is unbounded in 2tt while 
in the later case the phase value is fixed (Fig.QJ. Thus we 
may conclude that quantum synchronization gives a sig- 
nificant reduction of quantum fluctuations. At the same 
time fluctuations in p characterized by a p are not sensitive 
to the choice of /. In the limit of small h the numerical 
data of Fig|2|clearly show that the wave packet dispersion 
scales as _ 

VK. (5) 

In the synchronization regime the dispersion is even close 
to its minimum value with a x a p = az = h/2 (dashed 
line in Figl^J). At smaller dissipation 7 the asymptotic de- 
pendence a cx \f% starts at smaller values of ft (cf. left 
and right panels of Fig|3Jl . Indeed, the wave packet width 
grows dispersively during time 1 /j that makes a larger at 
small 7. The dependence © is in agreement with the data 
obtained in the regime of wave packet collapse in [231 f° r 
a smaller range of ft variation (note also [2*%]). 
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Fig. 4. (Color online) Poincare section for K = 0.8, 7 = 0.25 
at / = 1.1 (left) and / = 1.9 (right). Top: classical case; bot- 
tom: quantum case at ft = 0.05, points mark the average x,p 
positions of wave packet; light blue (gray) curves on x,p axes 
show at t — 5000 the quantum probability distribution in x 
and p respectively (arbitrary units) . Number of map iterations 
t = 5000, initial state is as in Fig. 



To check more accurately the scaling a x ,a p oc \fh we 
fix two values of force / and vary ft by three orders of mag- 
nitude (see FigEJ). One value of / is taken on a synchro- 
nization plateau when a classical attractor is a fixed point 



Other properties of quantum dissipative dynamics can 
be understood from Poincare phase space plots. In the 
classical case two typical examples are shown in FigQ] 
(top) with invariant curve and fixed point attractors. In 
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the corresponding quantum case we plot in the phase space 
the values of x, p averaged over a wave function of a given 
quantum trajectory at each map iteration (integer time t) 
(Fig0]bottom, quantum probability distributions are also 
shown). The quantum Poincare plot reproduces correctly 
the global structure of classical phase space but additional 
noise produced by quantum fluctuations is clearly visible. 
This quantum noise creates an effective width of an in- 
variant curve and replaces a fixed point by a spot of finite 
size. The size of these fluctuations is approximately given 
by a x and a p . They decrease with K according to the re- 
lation ©. Due to this quantum noise two quantum tra- 
jectories with the same initial state can drop on different 
attractors and thus contribute to different plateaus in the 
devil's staircase (Fig^). We indeed observed such cases 
simulating different quantum trajectories with the same 
initial state. However, convergence to different plateaus 
takes place only at relatively weak dissipation (7 = 0.05) 
when there are many different classical attractors close 
in a phase space and quantum noise may give transitions 
between them. At strong dissipation (7 = 0.25) the attrac- 
tors are well separated and different quantum trajectories 
converge to the same plateau. 

The above results show that the quantum synchroniza- 
tion has close links with the classical synchronization in 
presence of noise. The later has been intensively investi- 
gated and it is known that the synchronization plateaus 
are preserved if the noise amplitude is significantly smaller 
than their size [3] . In a similar way we may conclude that 
the quantum synchronization remains robust with respect 
to quantum fluctuations until their amplitude (propor- 
tional to <J x ,<7p oc y/h) remains smaller than the size of 
a synchronization plateau. We also checked that the clas- 
sical dynamics (j2Jl with additional noise in x,p, with dis- 
persion values <j x , o~ p taken from the quantum evolution, 
gives the dependence P vs. / which is close to the quantum 
result. This is illustrated in Fig. OD This correspondence 
remains valid if the size of the wave packet is sufficiently 
small (cr x ,a p <C 1). However, at large h or small 7 when 
o~x 1 a p ~ 1 this correspondence is broken and the quantum 
fluctuations cannot be reduced to effects of classical noise. 

Above we discussed quantum synchronization with ex- 
ternal periodic driving. It would be interesting to investi- 
gate the synchronization of two quantum nonlinear pendu- 
lum clocks to have a quantum version of the Huygens ex- 
periment 1 . Recently, first numerical simulations in such 
kind of regime (two quantum coupled Duffing oscillators) 
has been performed in |29j . Their results show that in 
the case when the quasi-integrable dynamics of oscillators 
becomes synchronized (entrained) the von Neumann en- 
tropy of one oscillator drops significantly. This result is in 
a qualitative agreement with our data showing that the 
dispersion of quantum state drops significantly on syn- 
chronization plateaus (Fig. |2J 30 . Further investigations 
of quantum synchronization in coupled nonlinear systems 
are of great interest. For example, quantization may pro- 
duce nontrivial effects in the synchronization transition in 
a disordered J J series array discussed in |3T| . 
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Fig. 5. Dependence of the average momentum P on static 
force / at K = 0.8, 7 = 0.25 (compare with the data in Fig.0 
left column). Right: quantum case at ft = 0.05 (same data as 
in Fig. left column). Left: classical simulation of classical 
dynamics with 200 classical trajectories in presence of noise in 
x,p with the dispersion of noise corresponding to the values 
cr x and <7p taken from the quantum trajectory simulation at 
h = 0.05 (right). 



Modern technological progress allows to study the regime 
of quantum synchronization with small size JJs (e.g. sim- 
ilar to those of [HI) or with cold or BEC atoms (e.g. like 
in |18| L This should open new perspectives for synchro- 
nization of quantum objects with possible applications to 
quantum computations. 

This work was supported in part by the EC IST-FET 
projects EDIQIP, EuroSQIP and (for OVZ) by RAS Joint 
scientific program " Nonlinear dynamics and solitons" . 
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